Revisiting Static and Dynamic Spin Ice Correlations in Ho2Ti207 
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Elastic and inelastic neutron scattering studies have been carried out on the pyrochlore magnet 
Ho2Ti2 07. Measurements in zero applied magnetic field show that the disordered spin ice ground 
state of Ho2Ti2 07 is characterized by a pattern of rectangular diffuse elastic scattering within the 
[HHL] plane of reciprocal space, which closely resembles the zone boundary scattering seen in its 
sister compound Dy2Ti2 07. Well-defined peaks in the zone boundary scattering develop only within 
the spin ice ground state below 2 K. In contrast, the overall diffuse scattering pattern evolves 
on a much higher temperature scale of ^ 17 K. The diffuse scattering at small wavevectors below 
[001] is found to vanish on going to Q=0, an explicit signature of expectations for dipolar spin ice. 
Very high energy-resolution inelastic measurements reveal that the spin ice ground state below ^ 
2 K is also characterized by a transition from dynamic to static spin correlations on the time scale 
of 10~^ seconds. Measurements in a magnetic field applied along the [110] direction in zero- field 
cooled conditions show that the system can be broken up into orthogonal sets of polarized a chains 
along [110] and quasi-one-dimensional f3 chains along [110]. Three dimensional correlations between 
P chains are shown to be very sensitive to the precise alignment of the [110] externally applied 
magnetic field. 

PACS numbers: 75.25.+Z, 75.40.Gb, 78.70.Nx 



I. INTRODUCTION 



Ho2Ti207 is a member of the rare-earth titanate, cu- 
bic pyrochlore materials, which crystallize into the FdSm 
space group^. This family of materials has been of great 
interest as the magnetic rare earth sites, the A-sites in 
the composition A2B2O7, reside on a network of corner- 
sharing tetrahedra, the archetypical structure for geo- 
metrical frustration in three dimensions^. The magnetic 
properties of this family of materials are very diverse, 
as the size and anisotropy of the relevant moments, as 
well as the nature of the coupling between moments on 
different sites, varies dramatically across the rare-earth 
series. Particular interest has focused on Ho2Ti207 and 
Dy2Ti2 07, both of which are characterized by strongly 
anisotropic, Ising magnetic moments with local [111] 
anisotropy, and net nearest-neighbor interactions which 
are ferromagnetic. The nature of the crystal field split- 
tings associated with Ho^+ in the Ho2Ti207 environment 
is such that large magnetic moments (~ 10 /i^ per Ho) 
are tightly constrained to point directly into or out-of 
the tetrahedra on which they reside^^^ 



This combination of a hard, local (111) anisotropy 
and coupled classical spins on the pyrochlore lattice has 
been well studied. Somewhat counterintuitively, anti- 
ferromagnetic near-neighbor interactions lead to a four- 
sublattice Neel state, which is not frustrated. However, 
an effective ferromagnetic near-neighbor interaction with 



local (111) anisotropy leads to a spin ice ground state, 
in which the spins on each tetrahedra are constrained 
by "ice rules" which require two spins pointing in and 
two spins pointing out on each tetrahedron^^^. There 
are six such degenerate configurations per tetrahedron, 
and hence an extensive number of configurations for an 
extended pyrochlore lattice. As is well known, this de- 
generacy is the same as that which arises from proton 
disorder in solid water ice, and accordingly this ground 
state is referred to as spin ice^. A unit cell of the py- 
rochlore lattice, containing 16 Ho^+ ions at the A-sites, 
is shown in Figure 1. The spins in Figure 1 have been 
chosen to obey the two-in, two-out ice rules for each 
tetrahedron, hence the spin configuration depicted is one 
of many degenerate spin ice ground states. 

The pyrochlore lattice can be decomposed in several 
fashions for convenience. For example, it can be thought 
of as an alternating stacking of triangular (edge-sharing 
triangles) and kagome (corner-sharing triangles) layers 
along the [111] direction^iSiii^iii. It can also be decom- 
posed into two orthogonal sets of chains, referred to as 
a and P chains, which run along the orthogonal [110] 
and [110] directions^ii^ii^ii^ii^. In this scenario, half of 
the spins in the system reside on the a chains, while half 
reside on the P chains. An illustration of this conceptual 
picture can be found in Figure 1, and as will become 
clear later on, such a decomposition can be very useful 
in the presence of a [110] magnetic field which breaks the 
crystal symmetry along one particular [110] direction. 
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FIG. 1: (Color online) Magnetic moments which satisfy the 
"ice rules" decorate a conventional unit cell of the pyrochlore 
lattice. The structure can be thought of as comprised of two 
orthogonal sets of chains running along the [110] and [110] 
directions. These are shown as red (dark) spins on the a 
chains and blue (light) spins on the orthogonal (3 chains. 



A useful starting point description for Ho2Ti207 is 
the following Hamiltonian for spins interacting on the 
pyrochlore lattice, known as the standard dipolar spin 
ice modeliii^ii^iii: 



(1) 



I > J 



Hj I ^ 



Hence the relevant ground state for Ho2Ti207 is ex- 
pected to be the spin ice state, a prediction which has 
now been well verified^^^^^^. 

Ho2Ti207 displays two low temperature anomalies in 
its heat capacity: a very low temperature feature which 
originates from nuclear hyperfine contributions, and a 
Schottky-like anomaly at ~ 2 K, which signifies the en- 
try to the spin ice state^^. The characteristics of the 
spin ice ground state in Ho2Ti207, as well as the analo- 
gous state in Dy2Ti207, have been of great interest for 
more than a decade. Recent theoretical work has focused 
on the possibility of emergent cluster correlations^ ^ai^, 
as well as the existence of magnetic monopoles^^ which 
fractionalize from the dipole moments originating from 
excitations out of the spin ice ground state. Recent ex- 
perimental work has investigated plateaus in the mag- 
netization of Ho2Ti207 and Dy2Ti207^^, and sought to 
identify the nature of the relevant field-induced phases 
and phase hounddiTie^^^^^^^^^ . A number of 
experiments have also probed the unusual dynamics of 
spin freezing in these materials^^-i2ii2^i2^i2ii2S. 

In this paper, we investigate several distinct themes. 
First we discuss the details of S(Q) as measured with 
time-of-flight neutron spectroscopy at low temperatures 
and compare these results to previous experimental 
and theoretical work. We focus specifically on topics 
not previously discussed: zone boundary scattering in 
Ho2Ti207 and the characteristic signatures of dipolar 
spin ice. We also investigate how S(Q) evolves with tem- 
perature as the spin ice state is destroyed upon warm- 
ing to T~ 2 K and above. We explore the dynami- 
cal behaviour of the Ho^+ spins as the spin ice state 
is approached upon cooling, and show direct evidence 
for spin freezing on a 10~^ second time scale. Finally, 
we characterize the decomposition of the spin ice state in 
Ho2Ti207 into polarized a chains coexisting with quasi- 
one-dimensional P chains in the presence of a [110] mag- 
netic field. 



The first term in Eq. 1 describes the relatively weak 
nearest-neighbor exchange interactions between Ho^+ 
moments, while the second term accounts for the Zee- 
man energy which arises due to the interactions of Ho^+ 
moments with the applied magnetic field. The third term 
describes the long-range dipolar interactions, which can 
be of considerable importance since the size of the mo- 
ments is very large. The S^^ notation indicates that the 
local Ising z-axis varies from site to site. Good estimates 
have been made for both J and D in Ho2Ti207, yielding 
values of -1.65 K (antiferromagnetic) and D~ 1.41 
K— As described by Bramwell and Gingras, the 
antiferromagnetic near- neighbor exchange in Ho2Ti207 
is sufficiently weak that it is overcome by the near- 
neighbor dipolar interactions, resulting in an effective 
near-neighbor interaction which is ferromagnetio^ai^ii^. 



II. EXPERIMENTAL DETAILS 

Several large single crystal samples of Ho2Ti2 07 were 
grown by floating zone image furnace techniques at Mc- 
Master University. The details of these growths were 
very similar to those which have previously been re- 
ported for Tb2Ti2 07^. The high quality single crystals 
of Ho2Ti207 which resulted from these growths were ~ 
5 grams in mass and had approximate cylindrical dimen- 
sions of 25 mm in length and 6 mm in diameter. A se- 
ries of neutron scattering measurements, utilizing both 
time-of-flight and backscattering techniques, were per- 
formed on two of these Ho2Ti207 single crystals at the 
National Institute of Standards and Technology Center 
for Neutron Research. In both of these experiments the 
samples were aligned such that the [HHL] plane in recip- 



3 



rocal space was coincident with the horizontal scattering 
plane. 

Time-of-flight neutron scattering measurements were 
performed using the Disk Chopper Spectrometer (DCS). 
DCS uses a series of seven choppers to create pulses 
of monochromatic neutrons, whose energy transfers on 
scattering are determined from their arrival times at the 
instrument's 913 detectors, located at scattering angles 
ranging from -30 to 140 degrees^^. Measurements were 
performed using both 5 A and 9 A incident neutrons. 
These measurements were used to map out the pattern 
of diffuse scattering in the [HHL] plane as a function of 
temperature and magnetic field, as well as to search for 
inelastic scattering with approximate energy resolutions 
of ~ 0.1 meV (5 A measurements) and 0.02 meV (9 A 
measurements). 

Much higher energy-resolution measurements were 
also carried out using the High Flux Backscattering 
Spectrometer (HFBS). In HFBS neutrons are backscat- 
tered from a spherically focusing monochromator, inter- 
act with the sample, and are backscattered a second time 
by a spherically focusing analyzer before reaching a bank 
of 16 detectors ranging from Q = 0.1 to 2.0 A-^ [Ref. 
31]. A Phase Space Transformation (PST) chopper is 
used to increase the neutron fiux of the instrument by 
Doppler shifting the incident neutron wavelength dis- 
tribution towards the appropriate backscattered wave- 
length. A high speed Doppler drive is used to oscil- 
late the monochromator, significantly expanding the dy- 
namic range of the instrument. In this experiment, in- 
elastic measurements were performed with a dynamic 
range of ± 17 /ieV and an approximate energy resolution 
of 0.92 jaeV. This is roughly a factor of 20 times greater 
energy resolution than the 9 A measurements collected 
with DCS. 



III. DIFFUSE MAGNETIC SCATTERING IN 
ZERO MAGNETIC FIELD 

Figures 2 a) and b) show S(Q) for Ho2Ti207 at T- 
0.2 K measured in the [HHL] plane. Figure 2 a) was 
collected using DCS and 5 A incident neutrons and 2 b) 
with DCS and 9 A incident neutrons, allowing smaller 
|Q| scattering to be accessed. The time-of-fiight neutron 
scattering technique simultaneously measures both the 
Q and energy dependence of the scattering. However, 
the diffuse scattering in both Figures 2 a) and b) is found 
to be resolution-limited in energy, and hence is elastic in 
nature. The diffuse scattering can be seen to peak up 
into well-defined rectangular regions, the most promi- 
nent of which are centered on the (0,0,1) and (0,0,3) 
positions in reciprocal space. The main region in the 
field of view which is devoid of diffuse scattering is cen- 
tered on (0,0,2) and can be seen to have an approxi- 
mately hexagonal shape. The overall pattern closely re- 
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FIG. 2: (Color online) Elastic scattering maps of S(Q) are 
shown for Ho2Ti2 07 at T=0.2 K. Measurements in a) em- 
ployed 5 A incident neutrons and elastic scattering cor- 
ersponds to AE<0.1 meV. b) employed 9 A incident neutrons 
and elastic scattering corresponds to AE<0.02 meV. The 
"x"s indicate the Q positions at which high energy-resolution 
backscattering measurements were performed. Dashed white 
lines indicate the dark angle caused by the sample environ- 
ment. 



sembles the zone boundary scattering which has been 
reported for Dy2Ti207^^, suggesting that such scatter- 
ing may be characteristic of spin ice as manifested in 
both of these materials. In Dy2Ti207, this zone bound- 
ary scattering has been attributed to either emergent 
hexagonal clusters^-, or to expectations from a gener- 
alized dipolar spin ice model which includes further- 
neighbour exchange interaction^. The measured S(Q) 
for Ho2Ti2 07 in Figure 2 differs in detail from both 
earlier reported neutron diffraction measurements on 
Ho2Ti207^^, and from the expectations of the near- 
neighbor spin ice model and the simple dipolar spin ice 
model (Eq. 

One of the most striking differences between calcu- 
lations for the near-neighbor spin ice model (with net 
near-neighbor ferromagnetic exchange and hard (111) 
anisotropy) and the simple dipolar spin ice model (which 
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includes the full Hamiltonian shown in Eq. 1) is a pro- 
nounced depression of the diffuse magnetic scattering 
near the origin of reciprocal space, at Q=0^^^^. This 
region of negligible scattering intensity at small Q arises 
due to the presence of dipolar interactions, which pre- 
clude the formation of a state with non-zero magnetiza- 
tion in zero applied field. This feature is unmistakably 
visible in our new measurements on Ho2Ti2 07, shown 
in Figure 2 b), which clearly reveal that the diffuse 
scattering drops off to background for L<0.6 for Q= 
[HHL]. The direct observation of this low Q depression 
in scattering intensity provides compelling evidence for 
the dipolar nature of the correlations in the ground state 
of Ho2Ti2 07. This feature of S(Q) can only be observed 
in the present set of measurements because they employ 
long wavelength 9 A incident neutrons. 

Zone boundary inelastic scattering has been observed 
in the spinel ZnCr204— where it has been interpreted 
in terms of strongly correlated spins in approximately 
independent hexagonal clusters which reside within the 
kagome planes of the pyrochlore structure. Similar scat- 
tering features have also been identified in the spinels 
CdFe2 04^^ and CdCr204^^. Theoretical work on a gen- 
eralized dipolar spin ice model has shown that elastic 
zone boundary scattering can also be accounted for by 
tuning exchange interactions, not included in Eq. 1, 
beyond near-neighbors^^. The importance of further- 
neighbor exchange interactions in spin ice has also been 
demonstrated by studies in which up to third nearest- 
neighbour exchange couplings have been used to explain 
the critical temperatures of observed field-induced tran- 
sitions in Dy2Ti2 07^'^'^. In any case, it seems clear 
that the nature of the zone boundary scattering is char- 
acteristic of a host of such highly frustrated pyrochlore 
magnets. 

S(Q) measurements using DCS with 5 A incident neu- 
trons, similar to those shown in Figure 2 a), were per- 
formed as a function of temperature up to T=20 K in 
order to investigate the evolution of spatial correlations 
across the spin ice transition near 2 K. As mentioned 
earlier, the entrance to the spin ice ground state is iden- 
tified by a Schottky-like anomaly in the heat capacity, 
so no conventional phase transition is expected. An 
additional surprisingly high temperature scale has also 
been identified - at least in Dy2Ti207. AC suscepti- 
bility measurements on Dy2Ti207 reveal a frequency- 
dependent anomaly in the 15 K regime^^^^^. A corre- 
sponding anomaly has also been identified in Ho2Ti207, 
but only in the presence of a DC field^^. This relaxation 
process is believed to be a thermally activated single-ion 
process, with an activation energy which is governed by 
the gap between the ground state and the first excited 
crystal field state^*^. The observed spin fluctuation 
rate is also typical for an individual spin flip process. It 
is the presence of these anomalies in the AC susceptibil- 
ity that determined the temperature scale on which we 



chose to explore the evolution of S(Q) in Ho2Ti207. 

Representative S(Q) measurements for IIo2Ti207 in 
the [HHL] plane at T=2 K (close to the spin ice Schottky- 
anomaly), T=4 K (well above the spin ice anomaly), and 
T=20 K (above the high temperature scale identified in 
the AC susceptibility measurements) are shown in Fig- 
ures 3 a), b), and c), respectively. These maps can be 
compared to the T=0.2 K data provided in Figure 2 a) to 
show how the zero field S(Q) in Ho2Ti207 develops with 
temperature, from an order of magnitude below the spe- 
cific heat anomaly at T ~ 2 K to an order of magnitude 
above it. As can be seen from these maps, the evolution 
of the measured S(Q) is somewhat subtle. This is not 
so surprising as the disordered spin ice ground state at 
T=0.2 K already displays exceptionally short correlation 
lengths (with spins correlated over a single tetrahedron), 
and hence warming up into the conventional paramag- 
netic phase cannot substantially shorten the overall spin 
correlations. The integrated intensity of line scans taken 
through the peak in the diffuse scattering at (0,0,3) is 
shown in Figure 3 d). This shows that the integrated 
diffuse scattering falls off gradually with slight upwards 
curvature on a relatively high temperature scale of ~ 
17 K - very close to the 15 K anomaly identified in the 
AC susceptibility measurements. This sets the scale on 
which S(Q) can evolve, since beyond ~ 17 K there is 
very little structure to S(Q) within our field of view. It 
should be noted that neutron spin echo measurements 
have detected evidence of diffuse magnetic scattering in 
Ho2Ti207 up to temperatures as high 800 K, although 
any indications of spatial correlations appear to vanish 
above --55 K^^. 

A collection of representative line scans taken through 
the zone boundary scattering are provided in Figure 
4. These scans have been taken along the [-2H+4/9, 
-2H+4/9, H+16/9] direction, allowing them to cut 
through the zone boundary scattering at its sharpest 
points, near (2/3, 2/3, 5/3) and (-2/3, -2/3, 7/3). The 
orientation of these cuts with respect to the overall pat- 
tern of diffuse scattering is displayed in the inset to Fig- 
ure 4 a), where the direction of the line scans can be 
identified by the dashed magenta line from (1,1,3/2) to 
(-1,-1,5/2). As a guide-to-the-eye, the hexagonal shape 
of the zone boundaries has also been highlighted by a 
series of dashed white lines. The line scans provided in 
the main panel of Figure 4 a) illustrate how the zone 
boundary scattering evolves with temperature in zero 
field, with data sets collected at T=0.2 K, 2 K, and 20 K. 
These cuts clearly demonstrate that well-defined peaks 
in the diffuse zone boundary scattering can only be iden- 
tified below the spin ice anomaly at T~ 2 K. Figure 4 
b) shows similar scans collected in the spin ice ground 
state of Ho2Ti207 at T=0.2 K, with a magnetic field 
applied along the [110] direction. As will be discussed 
in the following section, the application of a [110] field 
decouples the system into orthogonal a and (3 chains. 



5 




FIG. 3: (Color online) S(Q) for Ho2Ti207 in the [HHL] plane 
is shown at zero magnetic field for a) T=2 K, b) T=4 K, and 
c) T=20 K. The temperature evolution of the diffuse scat- 
tering at (0,0,3) is shown in d). Diffuse scattering intensities 
have been obtained by integrating over cuts taken along the 
[HH3] and [OOL] directions. The dashed line provided in d) 
is intended as a guide-to-the-eye, illustrating the plateau in 
diffuse scattering intensity above ~ 17 K. The error bars in d) 
and in all subsequent figures represent the standard deviation 
in the measurement. 



destroying the spin ice ground state at sufficiently high 
field strengths. 

An interesting subject which has recently attracted 
considerable attention is the potential presence of "pinch 
point" scattering in Ho2Ti207. Pinch point features 
have been theoretically predicted to arise in frustrated 
magnets as a consequence of enforcing the local ice 
rules^^^^^'^. These features, which are relatively sharp 
in specific directions, could be expected to broaden sig- 
nificantly upon passing through the spin ice anomaly. 
Sharp pinch point features in S(Q) have already been 
identified in Ho2Ti207 in the presence of a [111] mag- 
netic field^^. Such an applied magnetic field gives rise 
to a magnetization plateau which signifies the stabil- 
ity of the kagome ice phase of Ho2Ti207^^^^. As pre- 
viously mentioned, the pyrochlore lattice can be decom- 
posed into a series of alternating triangular and kagome 
planes which are stacked normal to the [111] direction. 
The kagome ice phase is characterised by spins which 
are locally constrained on the kagome planes, and spins 
which are polarized by the applied field on the triangular 
planes. 

The pinch point scattering predicted by the dipolar 
spin ice model is expected to arise at Brillouin zone cen- 
ters, such as (0,0,2). However, as the map of S(Q) pro- 
vided in Figure 2 a) clearly indicates, the (0,0,2) posi- 
tion is situated at the center of a region which is almost 
entirely devoid of diffuse scattering. Thus the present 
set of measurements reveal no indications of pinch point 
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FIG. 4: (Color online) Representative cuts taken through the 
sharpest points of the zone boundary scattering, at (2/3, 2/3, 
5/3) and (-2/3, -2/3, 7/3), showing a) the temperature evolu- 
tion of the scattering at /ioH=0 T, and b) the magnetic field 
evolution of the scattering at base temperature (T ^ 0.2 K). 
All cuts are taken along the [-2H+4/9, -2H+4/9, H+16/9] 
direction. The orientation of these cuts with respect to the 
[HHL] plane is illustrated by the dashed magenta line in the 
inset to panel a) . The zone boundaries have been outlined in 
white as a guide-to-the-eye. 



correlations in Ho2Ti207 at zero field, a result which 
is fully consistent with the previous work of Fennell et 
al^^. Unfortunately, due to the presence of a large field- 
dependent Bragg peak at (0,0,2) and the emergence of 
rod-like scattering along [OOL], we are unable to offer 
any definitive comment on the potential for pinch point 
features in an applied [110] field. Any evidence for the 
relatively subtle pinch point correlations is likely to be 
obscured by these much stronger field-dependent scat- 
tering features. 

It is also interesting to investigate how signatures of 
the entrance into the spin ice ground state are mani- 
fested in the spin dynamics of Ho2Ti2 07. As discussed 
earlier, below 2 K the energy-dependence of the diffuse 
magnetic scattering in Ho2Ti2 07 appears to be entirely 
elastic on the energy resolution scale which is accessible 
with DCS (0.1 meV with 5 A incident neutrons; 0.02 
meV with 9 A incident neutrons). Interestingly, this is 
very different from the diffuse magnetic scattering char- 
acteristic of the spin liquid ground state in Tb2Ti2 07 in 
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FIG. 5: (Color online) a) Representative energy cuts taken 
through zero field data sets at T=0.2 K and 20 K. Each cut 
has been binned over H= [-0.4,0.4] and L= [2.6,3.4] in order to 
capture the diffuse scattering at (0,0,3). b) The temperature 
evolution of the total elastic and inelastic scattering intensity, 
obtained by integrating over energy scans of the form shown 
in panel a). 



zero field^ i^^i^^ , which is inelastic on the energy scale of 
~ 0.3 meV^^. As shown in Figure 5 a), the energy width 
of the diffuse scattering can be seen to broaden slightly 
above 2K, providing an indication that spins are becom- 
ing more dynamic outside of the spin ice state. This is 
confirmed by the data in Figure 5 b), which shows that 
the gradual decrease in the intensity of the elastic diffuse 
scattering at higher temperatures is accompanied by a 
simultaneous increase in the intensity of the inelastic dif- 
fuse scattering. This shows that the static spins of the 
spin ice state undergo slow fluctuations as the temper- 
ature of the system rises above 2K. Improving on this 
energy resolution requires either the use of backscatter- 
ing techniques to measure S{Q^hu) directly, or neutron 
spin echo (NSE) measurements of S(Q,t) which have al- 
ready been carried out on Ho2Ti207. NSE measure- 
ments on Ho2Ti207 show a continuous evolution of the 
time-correlations for spins within the time window to 
which NSE is sensitive, over a temperature scale from 
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FIG. 6: (Color online) Representative backscattering energy 
scans taken at Q= (0,0,1) in Ho2Ti2 07. The scans displayed 
here were collected at T=1.4 K (within the spin ice state) 
and T=5 K (well above the spin ice state). 



0.3 K to 150 K^^^^. 

We have performed backscattering measurements near 
two wavevectors in the [HHL] plane; at (0,0,1), which co- 
incides with a peak in the pattern of diffuse scattering, 
and at (2/3, 2/3, 1/3), which lies near one of the nar- 
rowest regions of zone boundary scattering. These two 
positions are indicated by the pair of "x"s provided in 
Figure 2 a). Characteristic backscattering energy scans 
for Ho2Ti207 at T=1.4 K and T=5 K are shown in Fig- 
ure 6. As previously observed for the DCS data in Figure 
5, these representative scans clearly indicate a broaden- 
ing of the energy width above the Schottky anomaly at 2 
K. In order to extract quantitative values for the energy 
width of the diffuse scattering these data sets were fit to 
Lorentzian lineshapes which were convoluted by a Gaus- 
sian resolution function. The resulting FWHM of the 
scattering, which is characteristic of an inverse lifetime 
r, is plotted as a function of temperature in Figure 7. 
We observe that the diffuse scattering only approaches 
our resolution limit on entering the spin ice ground state 
below ~ 2 K. The paramagnetic state above T=2 K is 
characterized by a finite spin relaxation rate which is 
roughly temperature independent from ~ 3 to 30 K and 
increases with temperature above 30 K. Consequently, 
the distinguishing characteristics of the disordered spin 
ice ground state below the Schottky anomaly at 2 K are 
that the structure of the zone boundary scattering within 
S(Q) is well developed, and the structure is static on a 
time scale of 10~^ seconds. Neither of these statements 
is true for Ho2Ti207 at temperatures of T=2 K or above 
in zero applied field. 

In the inset to Figure 7 we see the temperature evolu- 
tion of the characteristic spin relaxation time r plotted 
on a log-log scale. Over the larger temperature window 
provided in this inset we can see three distinct relax- 
ation regimes: (i) the high temperature Arrhenius-like 
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FIG. 7: (Color online) Parameters extracted from fits to 
backscattering data obtained at the (0,0,1) and (2/3, 2/3, 
1/3) positions in reciprocal space. Shown in the main panel is 
the temperature evolution of the FWHM of the quasi-elastic 
signal. Provided in the inset is the temperature dependence 
of the relaxation time, r, which is inversely proportional to 
the width in energy. The spin ice state in Ho2Ti207 is char- 
acterised by spins which are static on the 10 ~^ second scale. 



relaxation above 30K, (ii) the temperature independent 
so-called quantum relaxation from 3K to 30K, and (iii) 
the diverging relaxation times indicative of spin-freezing 
below 2K. Plotted in this fashion, our data above 2 K 
fall into excellent agreement with previously reported 
NSE measurements on Ho2Ti207^^. However, what is 
unique to the present Ho2Ti207 backscattering mea- 
surements is our observation of the dramatic increase 
in relaxation times which occurs below the specific heat 
anomaly at 2 K. This behavior is fully consistent with 
the observed temperature dependence of the characteris- 
tic spin relaxation time in Dy2Ti207 as measured by AC 
susceptibility"^^, with only a slight shift in the relevant 
temperature scales distinguishing the two materials. 

Our backscattering measurements can also provide in- 
formation about how the integrated intensity of the dif- 
fuse scattering evolves with temperature. In particu- 
lar, the backscattering data collected at the (0,0,1) posi- 
tion are a useful complement to the DCS measurements 
shown in Figure 3 d) which describe the temperature 
dependence of the scattering at (0,0,3). Both (0,0,1) 
and (0,0,3) correspond to peaks in the overall pattern of 
diffuse scattering, and the consistency between the two 
sets of measurements is very good. In both cases the 
integrated intensity of the diffuse scattering appears to 
drop off on a much higher temperature scale than the 
relaxation time, somewhere in the vicinity of 15 to 20 K. 



IV. SCATTERING UNDER H||[110] 
ZERO-FIELD-COOLED CONDITIONS 

It has been appreciated for some time that the ap- 
plication of a magnetic field along [110] and equivalent 
directions should decouple the spin ice ground state into 
polarized a chains (which point along the field) and per- 
pendicular P chains^ii^ii^ii^ii^, as illustrated in Figure 
1. This decoupling has been demonstrated directly in 
Dy2Ti207 by neutron scattering experiments^^. Related 
measurements have also been performed on Ho2Ti207— , 
although to date maps of the low temperature diffuse 
scattering have only been reported for Dy2Ti2 07. 

In Figure 8 the measured S(Q) for Ho2Ti207 in the 
[HHL] plane is shown at base temperature (T=0.2 K) 
for H II [110]. The data sets provided in these four pan- 
els correspond to applied magnetic field strengths of a) 
0.2 T, b) 0.4 T, c) 0.9 T, and d) 2.5 T. When com- 
bined with the zero field S(Q) data displayed in Figure 2, 
these maps illustrate the evolution of the diffuse scatter- 
ing from broad, rectangular, zone boundary scattering at 
zero field to narrow, rod-like scattering extending along 
[OOL] at finite fields. This rod-like scattering can be in- 
terpreted as arising from sheets of quasi-one-dimensional 
scattering from decoupled [3 chains intersecting with the 
[HHL] scattering plane. Coinciding with the develop- 
ment of this rod-like scattering is the growth of the scat- 
tering intensity at the (0,0,2) Bragg peak, which is in- 
dicative of the polarization of the a chains. Although it 
is possible to have ordered magnetic states with polar- 
ized a chains and negligible (0,0,2) peak intensity, in the 
case where there is no net moment within the (5 chain 
sublattice, the intensity of the (0,0,2) peak can be taken 
as a direct measure of the a chain polarization. For 
these measurements, and all subsequent measurements 
reported below, the magnetic field was applied at low 
temperatures following a zero-field cooled protocol. 

Cuts through the diffuse scattering shown in Figure 8 
were taken along the [OOL] direction, parallel to the rod 
of scattering, and along the [HH3] direction, perpendic- 
ular to the rod, as shown in Figures 9 a) and b). These 
cuts were then fit to an anisotropic Lorentzian using an 
Ornstein-Zernike form, for the purpose of extracting cor- 
relation lengths parallel and perpendicular to the rods. 
The resulting correlation lengths in real space are shown 
as a function of field in Figure 9 d). These fits show that 
the diffuse scattering is fully three-dimensional in zero 
field, as expected, with identical correlation lengths of ~ 
3.6 A, or one near-neighbor separation, in all directions. 
The scattering immediately becomes anisotropic upon 
the application of a [110] magnetic field, with both sets 
of correlation lengths growing as a function of increas- 
ing field strength. As can be seen in Figure 9 d), the 
transition to a quasi-one-dimensional structure is com- 
plete by ~ /ioH=l T, at which point the correlations 
along the [3 chains are resolution-limited (> 100 A) and 
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FIG. 8: (Color online) S(Q) for HoaTisOr in the [HHL] plane 
is shown for T=0.2 K and H|| [110]. The maps displayed above 
correspond to field strengths of a) 0.2 T, b) 0.4 T, c) 0.9 T 
and d) 2.5 T. In each case the magnetic field was applied 
following a zero-field cooled protocol. 



those between /3 chains are saturated at 10 A, or roughly 
two interchain distances. The integrated intensity of the 
(0,0,2) Bragg scattering was also monitored as a func- 
tion of [110] magnetic field, as shown in Figure 9 c). As 
noted earlier, the (0,0,2) intensity can be interpreted as 
a measure of the polarization of spins on the a chains, 
and it is observed to grow sharply and approach satura- 
tion by 0.9 T. This corresponds very closely with the full 
development of the quasi-one-dimensional correlations of 
spins residing on the P chains, a result which agrees very 
well with previous neutron results^^. 

We now return to the zone boundary scattering shown 
in the line scans of Figure 4, and specifically how these 
spin ice correlations evolve with the application of a 
[110] magnetic field. The field dependence of the zone 
boundary scattering is shown in Figure 4 b) for T=0.2 
K and representative field strengths of 0, 0.1, 0.3, 0.5, 
and 0.9 T, respectively. These scans capture the peaks 
from the narrowest regions of the diffuse scattering, near 
(2/3,2/3,5/3) and (-2/3,-2/3,7/3), as well as the field- 
dependent Bragg scattering observed at (0,0,2). We find 
that the intensity of the zone boundary scattering begins 
to drop with the application of any finite [110] magnetic 
field, with evidence of distinct peaks vanishing by ~ 0.5 
T. This coincides with the point at which the (0,0,2) 
Bragg scattering, representing the polarization of spins 
on the a chains, is almost fully developed. The complete 
elimination of the diffuse scattering along the cuts shown 
in Figure 4 b) requires a field of 0.9 T, at which point 
the intensity of the (0,0,2) peak is completely saturated. 

We also note that the scattering intensity of the (0,0,2) 
Bragg reflection is weak, but unmistakably non-zero, 
even in the absence of an applied magnetic field. This 



FIG. 9: (Color online) Representative cuts of the data shown 
in Fig. 8 taken along a) [OOL] and b) [HH3] are shown at 
T=0.2 K for applied magnetic fields of /ioH=0, 0.5, and 2.5 
T. The field dependence of the (0,0,2) Bragg peak intensity is 
shown in c). The correlation lengths of the diffuse scattering 
along [OOL] and [HH3] are shown in d), where the dashed line 
indicates the resolution limit of the instrument. The quasi- 
one- dimensional correlations of the f3 chains are established 
by a - 1 T magnetic field ||[110]. 



(0,0,2) scattering can clearly be observed in the zero-field 
measurements of S(Q) provided in Figures 2, 3, and 4. 
This weak Bragg feature is structurally forbidden within 
the FdSm space group appropriate to the pyrochlores, 
and appears to be temperature independent up to at 
least T=20 K. Earlier neutron diffraction measurements 
on pyrochlores such as Tb2Ti207^^'^^ and Er2Ti2 07^^ 
have identified similar scattering at the (0,0,2) position, 
but in general there has been a tendency to associate 
this feature with harmonic Bragg (A/n) contamination. 
Since the present time-of-flight measurements (like those 
of Refs. 39 and 44) do not employ Bragg reflection from 
a monochromator cystal, we can be assured that such 
scattering is not due to higher order contamination, and 
likely indicates a weak departure from the perfect py- 
rochlore structure in Ho2Ti207 and perhaps in all real 
pyrochlore materials^. 

A series of S(Q) measurements were also performed 
for Ho2Ti207 with the applied magnetic field offset from 
[110] by a rotation of 7.6 ib 3 degrees about the [OOL] di- 
rection. These measurements should be contrasted with 
the data discussed above for which the magnetic field is 
nominally applied precisely along [110], but in practice is 
aligned to within an accuracy of ± 1 degree. The identi- 
fication of the precise misalignment of the magnetic field 
comes about by virtue of the relative positions at which 
several Bragg peaks were observed on the three detector 
banks of DCS. In this case, the appearance of the (0,0,2) 
Bragg reflection on the middle detector bank, combined 



with allowed FdSm Bragg peaks such as (1,1,3) and (-1,- 
1,3) appearing on the upper and lower detector banks, 
is consistent with a field misaligned from [lIO] by 7.6 
± 3 degrees about [OOL]. Note that the same rotation 
which changes the orientation of the applied field will 
also affect the scattering plane used in the experiment. 
If the applied field direction is offset from [110] by some 
angle 6 about [OOL], then the new scattering plane can 
be expressed as: 

[H{cosO + sinO) , H{cosO - sinO) , L] (2) 

Thus, for a field offset of ^ ~ 7.6 degrees, the scatter- 
ing plane should be roughly coincident with the [1.13H, 
0.86H, L] plane in reciprocal space. 

A comparison of the measured S(Q) for Ho2Ti207 
in the presence of aligned and misaligned applied mag- 
netic fields is provided in Figure 10. Both of the data 
sets depicted here were collected using the same crys- 
tal of Ho2Ti2 07, and both sets of measurements were 
taken at T=0.2 K and /ioH=0.5 T. Figure 10 a) shows 
the nominally aligned data, where we can see quasi- 
one-dimensional rods of scattering along [OOL] , but clear 
short-range correlations along the [HHO] /3 chain di- 
rection, consistent with the data shown in Figure 4. 
In addition, one can clearly observe the strong (0,0,2) 
Bragg scattering which arises due to the polarization of 
a chains. Figure 10 b) shows the same data set with the 
applied field canted by 7.6 ± 3 degrees with respect to 
the [OOL] direction. Cuts along the [OOL] direction for the 
data sets shown in Figures 10 a) and b) are overlaid in 
c). It is evident from b) and c) that the diffuse scattering 
along [OOL] is far more one-dimensional in the presence of 
the canted magnetic field, with much less structure vis- 
ible in both the S(Q) maps and the cuts taken through 
the data. Interestingly, while the properties of the [OOL] 
diffuse scattering are significantly altered by field mis- 
alignment, the nature of the Bragg scattering at (0,0,2) 
remains essentially unchanged. 

The effect of the canted [110] field can be understood 
in terms of its decomposition into a large aligned [lIO] 
field along the a chains, and a smaller field perpendic- 
ular to the a chains due to misalignment. In this case, 
the perpendicular component of the field is aligned along 
[110], parallel to the (3 chains. The weak, primarily dipo- 
lar, interactions between /3 chains are expected to be 
antiferromagnetic in nature, and this is consistent with 
the pattern of diffuse scattering observed in Figure 8. 
The application of a small uniform magnetic field along 
the direction of the /3 chains would frustrate these an- 
tiferromagnetic correlations, resulting in the enhanced 
one-dimensionality which we observe in the presence of 
the canted field. We also expect that small offsets to the 
field direction should have little measureable effect on 
the polarization of the a chains, since the component of 
the field along [110] will be only slightly diminished by 
misalignment. Once again, this prediction appears to be 
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FIG. 10: (Color online) A comparison is made between diffuse 
scattering with a) H applied || [110], and b) H canted with 
respect to [110] by 7.6 ± 3 degrees. The canted field implies 
a field component along the f3 chains which frustrates their 
short range ordering and results in a more-one-dimensional 
structure. The vertical axes in a) and b) are a function of 

0, the angle between [110] and the direction of the applied 
field, which can be written as [H(cos^+sin^), H(cos^-sin^), 

01. The cuts provided in c) have been taken along the [OOL] 
direction, through the diffuse scattering shown in a) and b). 



fully consistent with the data shown in Figure 10. 

One implication of this result is that such offsets to 
the field orientation, even in the case of nominal "per- 
fect" alignment, may be sufficient to frustrate the three 
dimensional correlations between P chains. In turn, this 
frustration may cause the spin ice system to remain in a 
quasi-one-dimensional short-range ordered state, rather 
than arriving at its true equilibrium state^^, character- 
ized by long-range ordered P chains. The significance 
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of small deviations from perfect alignment has also been 
demonstrated theoretically by Monte Carlo studies on 
dipolar spin ice^^, where misalignments of only one de- 
gree have been shown to suppress antiferromagnetic in- 
terchain order in modest [110] applied fields. 

V. CONCLUSIONS 

In conclusion, we have carried out neutron scattering 
measurements on the spin ice ground state of single crys- 
tal Ho2Ti2 07 as a function of temperature and applied 
magnetic field. In zero magnetic field, S(Q) at T=0.2 
K is characterized by diffuse, elastic scattering which 
closely resembles the zone boundary elastic scattering 
seen in Ho2Ti207's sister spin ice material, Dy2Ti207. 
We have explicitly shown that the diffuse scattering falls 
to zero as |Q| goes to zero, consistent with theoretical 
predictions from the dipolar spin ice model. We have 
also shown that the diffuse zone boundary scattering de- 
velops a well-defined structure only within the spin ice 
state for T< 2 K. Furthermore, it has been shown that 
the integrated intensity of the strongest diffuse features, 
the rectangles centered at (0,0,1) and (0,0,3), evolves on 
a surprisingly high characteristic temperature scale of 
~ 17 K, similar to that of the 15 K anomaly observed 
in Dy2Ti207. Very high energy-resolution inelastic neu- 
tron scattering measurements have demonstrated that 
the spin ice ground state of IIo2Ti207 is static on a time 



scale of 10 ^ seconds. 

We have also measured the evolution of S(Q) at T=0.2 
K for IIo2Ti207 in the presence of a [110] magnetic 
field. Application of the field decomposes the system 
into spins residing on polarized a chains and those re- 
siding on weakly correlated /3 chains, leading to the for- 
mation of a quasi-one-dimensional magnetic substruc- 
ture. We have shown that the peaks in the diffuse scat- 
tering at the zone boundaries vanish for magnetic field 
strengths which substantially polarize the a chains. Fi- 
nally, we have demonstrated that the nature of the three- 
dimensional correlations between quasi-one-dimensional 
/3 chains is very sensitive to the precise alignment of the 
externally applied [110] magnetic field. We hope that 
this work will help to guide and inform future studies of 
the spin ice state in the rare earth pyrochlores. 
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